Point vortices on the sphere: a case with oppo- 
site vorticities 
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Abstract. We study systems formed of 2N point vortices on a sphere with N 
vortices of strength +1 and N vortices of strength —1. In this case, the Hamil- 
tonian is conserved by the symmetry which exchanges the positive vortices with 
the negative vortices. We prove the existence of some fixed and relative equihb- 
ria, and then study their stabiUty with the "Energy Momentum Method" . Most 
of the results obtained are nonlinear stability results. To end, some bifurcations 
are described. 
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1 Introduction 

The study of the motion of many eddies (vortices) started one century ago when 
Helmhotz Q| introduced the point vortex model. With progress in dynamical 
systems in recent decades, the iV-vortex problem has re - appeared. Numerous 
papers have been written on vortices in the plane |Ar82 , Ar83a , ArSSb , AV98 



LR9£, 3a92 |. From a geostrophic point of view, it is interesting to study point 



vortices on the sphere. In this paper we consider a non-rotating sphere, that is 
we focus only on the spherical geometry of the Earth. In a forthcoming paper we 
will consider a rotating sphere and then deal with the forced symmetry breaking 
theory |CLOO|. As many results in this paper depend only on the phase space 
and the symmetries of the Hamiltonian, the study of point vortices on the sphere 
may also be fruitful for other problems which arise in physics: charged particles 
on a sphere, a-Euler model in turbulence, motion of cells on a spherical flame, 
zeros of wavefunctions in quantum physics. . . 

The motion of three vortices of arbitrary strength on a sphere is studied in 
[ [KN98| and | |KN99| The stability of the relative equilibria given in | |KN9^ is 
computed in |PM9S] and numerical simulations are done in |MPS99 . Relative 



equilibria formed of N vortices are treated in |LMR01| and the stability of a 
latitudinal ring of N identical vorticities is computed in [LMR]: it appears that 
linear stability results obtained by [PD93| are in fact Lyapunov stability results. 

The dynamics of point vortices on a sphere is given by an Hamiltonian system 
(|B77|). The Hamiltonian is invariant under rotations of the sphere, reflections 
of the sphere and permutations of identical vortices. In this paper, we consider 
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A'' vortices of vorticity +1 and N vortices of vorticity —1, where one has a new 
symmetry for the Hamiltonian: the permutation which exchanges the positive 
vortices with the negative vortices. However, some of these symmetries (eg re- 
flections) are not symmetries of the equations of motion, they are time-reversing 
(or, anti-symplectic) symmetries (Section 2.2). From Noether's theorem, the ro- 
tational symmetry provides three conserved quantities, the components of the 
centre of vorticity map ^ :V ^M.^ where V is the phase space. 

In Section 3, we determine some equilibria and relative equilibria. Relative 
equilibria are orbits of the group action which are invariant under the flow and 
correspond of the motions of the point vortices which are stationary in a steadily 
rotating frame. In the same way in which equilibria are critical points of the 
Hamiltonian if, relative equilibria are critical points of the restrictions of H to 
the level sets 3'~^(/u). The fixed point set of a group of symmetries with some 
anti-symplectic symmetries might not be invariant under the flow contrary to 
those formed of symplectic symmetries. However, it follows from the Principle 
of Symmetric Criticality |P79, Mi71| that critical points of the restriction of 
H to such a fixed point set are critical points of the full Hamiltonian, and 
are therefore equilibrium points. An analogous statement holds for relative 
equilibria. Relative equilibria derived with these results do not depend on the 
expression for the Hamiltonian, they depend only on the symmetries, the phase 
space and the momentum map. This technique is especially useful for finding 
relative equilibria in high dimensional systems with large symmetry groups. It 
is probable that there are "asymmetric" relative equilibria but they can not be 
found by the symmetry adapted methods of this paper. See | AV98| for examples 
of asymmetric relative equilibria in the plane. 

In [LMR01|, the authors consider also anti-symplectic symmetries and use 
the Principle of Symmetric Criticality, but the approach is quite different. In 
the case of N identical vortices, the Hamiltonian goes to minus infinity when 
some vortices collide. Thus the Hamiltonian must have a maximum somewhere, 
so there exists an equilibrium poin t. Th ey generalized this idea to fixed point 
sets (see the remark after Theorem 3.10 for an example). 

Section 4 is devoted to analyzing the stability of the relative equilibria de- 
termined in Section 3. The suitable concept of stability for relative equilibria 
in Hamiltonian system is Lyapunov stability modulo a subgroup. The stability 
study is realized using the energy momentum method |Si9l| , |Ma92| , |Pa92| , |LS98| 



OR99| . This is stronger than linear stability. For configurations with a small 
number of vortices, there are some Lyapunov stable relative equilibria. For a 
large number of vortices the relative equilibria we study are all linearly unstable, 
and in general we do not know whether there are other stable relative equilib- 
ria. When a relative equilibrium is claimed to be linearly stable, this means in 
addition that the Lyapunov stability criterion (the energy momentum theorem, 
a kind of Dirichlet's criterion for relative equilibria) does not apply. Finally, we 
describe some pitchfork and Hamiltonian- Hopf bifurcations in systems of four 
and six vortices. 
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2 The point vortices system on the sphere 
2.1 Equations of motion 

Let N be an integer, N > 2. The phase space for the 27V-vortex problem is a 
product of N spheres where the diagonal has been removed in order to avoid 
collisions: 

V = {{xi,...,X2n) e X ■ ■ ■ X \ ^ Xj if i^j}. 

In the previous expression, each Xi represents a point vortex on the sphere, of 
vorticity A^. 

In principle we assume that half of the vortices are of vorticity +1 and 
the other half of vorticity —1, though this assumption can be relaxed when 
polar vortices are considered. We could write an element of V in two different 
manners: 

- either x = (xi, . . . ,xn,xn+i, ■ ■ ■ ,X2n) with Xj = +1 and X^+j = —1, j = 
l,...,iV, 

- or X — (xi , . . . , X]s[ ,x[,..., x'j^ ) with Xj = + 1 and = — 1 , j = 1 , . . . , iV. 

For notational convenience, arrows above vectors will be omitted most of the 
time. We call a f+Jvortice (resp. a (—)vortice) a vortice with vorticity +1 (resp. 

- 1), a (+)ring (resp. a (—)ring) a latitudinal regular gon of (+)vortices (resp. 
(— )vortices), and a (±)ring a latitudinal regular polygon formed of alternating 
(+)vortices and (— )vortices. Throughout this paper, Tq denotes the rotation of 
50(3) about the (Oz) axis with angle a, and Sx, Sy, Sz denote respectively the 
reflections of 0(3) of planes (yOz), (xOz), (xOy). 

The equations governing the motion of 2N vortices are given by the following 
system of 6A^ equations | B77 , KN9^ : 

x,=XHix).^J2 ^^^^^^ , z = l,...,2iV. 

.... -L Xj ■ Xj 

Since the vortices are constrained to lie on the sphere, these equations are not 
independent. 

Let 9i be the co-latitude and (pi be the longitude of Xi. In spherical coordi- 
nates, the system is formed of the following AN equations [ KN98{ : 

sm9^ sini (hj — (hi) 

2N 

■ n 1 o \ sin6'i cos6'j - sin6'j cos6'icos(0j - 0j ) . 
sm Oi (fn = 2 Xj — , z = 1, . . . , 2N 

where Ijj = 2(1 — cos6'i cosdj — s'mOi s'mOj cos((/)i — (j)j)) is the square of the Eu- 
clidian distance jjxi — Introducing the conjugate variables qi — \/\Xi\ cos9i 
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and Pi — sign{Xi)^J\\i\(f)i, and the Hamiltonian 

i/ = ^A,A,lnZ2. (2.1) 

we put the system in Hamihonian canonical form: 

dH 

dH 
OQi 

for all i = l,...,2N. 

The phase space has so a symplectic structure given by 

Ol! — Xi sin 9i dOi A d4>i 

i 

and a Poisson structure given by 

^ Xi \d cos 9i d(j)i d(t>i d cos 6i J 

For any function /, we have / ~ {f,H}. Embedding 5^ in M"^, the Poisson 
bracket becomes: 

where (a, 6, c) — {a x b) ■ c. 
2.2 Symmetries of the equations 



In this part, we apply the work of Lim, Montaldi, and Roberts [ LMROl ], Let r 
be the permutation of 5*2^, t — Y[i=i where Tij is the transposition which 
exchanges Xi and Xj. That is, r is a permutation of order two which exchanges 
(+) vortices with (—) vortices. 

Let G = 0(3) X Sn x Sn x 1.2[t] where H yi K is the semi-direct product 
of the groups H and K. The group G acts on V in the following manner: 

g-{xi, . . .,XN,x[, . . ■,x'^) = (Ax^fc^(i), . . . ,yla;^fc^(jv),^a;^fcCT'(i'), ■ • . , Ax^k„,(^ff,^) 

where g — {A, a,a' ,t''^) G G. Let x be the following homomorphism, 

X-G ^ {-!,+!} 
(A,(7,a',T'=) ^ {-!)'' dot A 

The following proposition gives symmetries of the Hamiltonian and the dynam- 
ical system. 
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Proposition 2.1 

• The Hamiltonian H is invariant under the action of G onV. 

• Elements ofKeT{x) 'J^g symmetries of Xh, while elements of G \ KeT{x) are 
time reversing symmetries of Xh ■ 

Proof, (i) Easy to check. 

(ii) Let g = {A, a, a' ,t^) gG and f,he G°^{V), we have 
Since {f,H}{x) = df^-XH{x), 

^nig ■ x) = xig) dgx ■ XH{x)n 



The action of G is semi-symplectic (see ||lVIoR | for the study of such an action 
in the reduced space): elements g g Ker(x) are symplectic, they preserve the 
symplectic form and the Poisson bracket, while elements of g G G \ Ker(x) are 
anti- symplectic, the symplectic form and the Poisson bracket are transformed 
into their opposite. The anti-symplectic elements are time reversing symmetries 
of Xh- if t I— s- x(t) is a solution of the dynamical system with x(0) = xq, then 
the solution with initial value g ■ xq is 1 g ■ x{—t). 

The isotropy subgroup Gx — {g G \ g ■ x = x} of a, point x ^ V will be 
called the symmetry group of x. Let K he a subgoup of G, we recall that the 
fixed point set of K is Fix(i^, T') ~ {x <^ V \ g- x = xyg^ K}, and that if 
K is formed of symplectic elements, then Fix(_fC, V) is invariant under the flow. 
When there are no ambiguities, we write Fix(iir) instead of Fix(ii', T'). 



Remark. We can construct a symplectic element from two anti-symplectic ele- 
ments (a reflection composed with the permutation r). For example, consider K 
the subgroup of G generated by the two elements {t^^, (1 . . . N)^^ , [!'... N')^^) 
and {sx,t). Fix(if) is the set of latitudinal semi-regular 2A^-gons formed of al- 
ternating vorticities. The subgroup K is formed of symplectic elements, Fix(i4r) 
is then invariant under the flow. Actually in this case the dynamic can be 
completely described: see [To]. 



2.3 Conserved quantities 

From Emmy Noether's Theorem, Hamiltonian systems with continuous symme- 
try group satisfy conservation laws. These conserved quantities are the compo- 
nents of the momentum map. 

Consider the diagonal action of 5*0(3) on V . Let ^ e so(3) and ^-p the vector 
field obtained by differentiating g-x with respect to g in the direction ai g — e. 
Identifing 5o(3) with equipped with the cross product as Lie bracket, we get: 

iv{x) = (^ X xi, . . . ,^ X a;2Ar) 
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Recall that a momentum map J : "P — > 5o(3)* ~ for this action is such that 

for all / e C°°{V) and ^ G so (3) ~ M^. We identify the natural bracket between 
so (3) and its dual with the canonical scalar product of R^. 
Let the center of vorticity vector be the following vector: 

2N 

Proposition 2.2 The momentum map for the action of S0(3) is J{x) = — $(a;). 

Proof. It follows from the definition of momentum maps J and the expression 
of ^-p that for all / € C°°{V), 

Thus for all C and j, Xj^ x xj = -^^^^^ x xj, therefore = -Xj for 

all J = 1, . . . , 2A''. Hence the momentum map is J{x) = — J2j ^j^j up to a 
constant. □ 

Consequently, the components of $ arc three conserved quantities (we could 
compute this directly differentiating $ with respect to time). 

Throughout this pciper, we will take ci frcinie (O, e^^ 6z ) of such that 
$ is parallel to the {Oz) axis. We obtain the following properties of $ from its 
formula. 

Proposition 2.3 

• ^ is G-equivariant with the following action on so(3)* ~ R^: 

g- ( = AC 

where g = {A, a, a', t^) G G. 

• Let K he a subgroup ofG. If x e Fix{K,V), then <^{x) G Fix{K,R^). 

3 Equilibria 

3.1 Fixed equilibria 

To emphasize the difference between equilibria and relative equilibria, we call 
the former fixed equilibria. Thus a fixed equilibrium is a critical point of the 
Hamiltonian. First, we give a general result which will serve to prove existence 
of fixed equilibria. 
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Theorem 3.1 Let K be a subgroup of G. An isolated point in Fix(_ftr) is a fixed 
equilibrium. 



This Theorem was given first by Michel [Mi71 , and it is also an application of 
the Principle of Symmetric Criticality of Palais [ |P79| ] which states that if the 
directional derivatives dHx{u) vanish for all directions u aX x tangent to Fix(iir), 
then directional derivatives in directions transverse to Fix(/'ir) also vanish. In 
particular, an isolated point of Fix(/'ir) is a critical point of H . 

Note that equilibria derived by this Theorem do not depend on the form 
of the Hamiltonian, they depend only on the symmetry group and the phase 
space. For example, the system of N point charges on a sphere is Hamiltonian 



with potential ( |GE92| ) 



where is the charge of the i-th particle and Uj is the Euclidian distance 
between particles i and j. Hence, fixed equilibria of the vortex system described 
thereafter are also fixed equilibria of the point charges system. But the stability 
results for fixed equilibria of the vortex system do not hold for the corresponding 
fixed equilibria of the point charges system, since the stability analysis depends 
on the expression for the Hamiltonian. 



Kidambi and Newton [KN98| found that a fixed equilibrium formed of vor- 
tices Xi with vorticities A; must statisfy J^ii^i ~ J2j = 0. In the present 
case of 27V vortices a^i, . . . ,X2n with Ai = . . . = \n = — Ajv+i = . . . = — A2Ar, 
this implies that the centre of the sphere is the isobarycentre of the configuration 
(bary centre with identical weights). 

Proposition 3.2 The following configurations are fixed equilibria: 

- an equatorial (±)ring. 

- a regular tetrahedron formed of 4 (+)vortices and its tetrahedron dual formed 
of 4 (—Jvortices. 



Following the notation introduced in |LMR01|, we denote the equatorial (±)ring 
by D2Nh{Re) (see figure 1). 

Proof. Let K be the symmetry goup of an equatorial (±)ring. We have tt{K) = 
D2N X '^2[sz] where tt is the cartesian projection n : G ^ 0{3) and Dpf is the 
dihedral group of order 2N. Consider also K' — {{g, (Jg^ , (Jg^), g & T) C G 
where T is the subgroup of 0(3) formed of the symmetries of a regular tetrahe- 
dron, and CTg is the permutation induced by g on the vertices. The reader can 
check that Fix(_ftr) and Fix{K') are formed of isolated points. So Fix(_ftr) and 



Fix(if') are formed of fixed equilibria by use of Theorem 3.1. □ 



Remark. This proposition is also true if the tetrahedron and its dual do not 
have opposite vorticities. Indeed, the proof does not use the permutation r 
which exchanges the (-I-) with the (— ). In the same manner, the following 
arrangements are fixed equilibria: a cube formed of 8 identical vortices and its 
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Figure 1: Equatorial (±)ring D2Nh{Re)- 



octahedron dual formed of 6 identical vortices, a dodecahedron formed of 20 
identical vortices and its icosahedron dual formed of 12 identical vortices. For a 
topological proof and some other fixed equilibria formed of regular polyhedra, 
see [|LMROlf. 



3.2 Relative equilibria 

Let Xe £ V and Xe(t) the dynamic orbit of Xh with a;e(0) — Xe- Then Xe 
is called a relative equilibrium if for all t there exists gt G S'0(3) such that 
Xe{i) = gt ■ Xe- 

Let n — ^{xe) G K'^- We recall that Xe is relative equilibrium if and only if 
there exists ^ G so(3) ~ M'^ such that critical point of 

H^{x)^H{x) + {^x)-^l,o 

where the pairing () between M'^ and its dual is the canonical scalar product of 
R^. We can say equivalently that relative equilibria are critical points of the 
restriction of H to <I>~^(/x) since the level set $~^(/^) are always non singular 
for point vortex systems of more than two vortices. 

The solution is given by Xe{t) = exp(t^) • Xe, the vortices turn therefore 
uniformly around ^ and ^ is called the angular velocity of the relative equilibrium 

Xe • 

We can choose the action of G on so (3) such that if a relative equilibrum 
has a particular symmetry, then its angular velocity has also this symmetry, as 
the following result shows. 

Proposition 3.3 Let Xe be a relative equilibrium with angular velocity f . Then 
^ satisfies 

where x is defined in Section \2.'J^ Moreover, ^ G Fix(Gj:^ , so(3)) where the 
action of G on so(3) is defined by g ■ ^ = xid) f . 
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Proof. Let a: be a relative equilibrium with angular velocity ^, and g E Gx- 
Hence, dHx — d^^)- Since dH^ = dHg.^ ■ dg^ = dH^ ■ dg^ and d^x • dgx = 
Xig) CoAdg d^x, 

{Adg C, d^x) - CoAdg-1 = ^(g){^^ d<S>x ■ dgx) = xi9){^, d<fx). 

Thus Adg^ = xig)^ since $ is a submersion. Clearly ^ g Fix{Gx,so{3)) with 
respect to the action g ■ £, — x{9) ^ 

The following property is due to the rotational symmetry. 

Proposition 3.4 The angular velocity of a relative equilibrium is parallel to 
the vortices rotate around the momentum vector. 

Proof. Let Xe be a relative equilibrium with angular velocity ^ and jj, = ^{xe). 
We have Xe{t) = exp(i^) • Xe € <i>~^(^) n G • Xe = • Xe, so exp(t^) e G^ and 
CoAdcxp(t5) IJ- — t^- As CoAdyi — A for A £ 50(3), we get exp(i^)/i — ^. A 
differentiation of the last identity will give the result. □ 



Remark. As pointed out by Kidambi and Newton [ KN98 1 , $ • depends only 
on the distances lij for all i, and so these are conserved quantities if xi, . . . , X2n 
is a relative equilibrium. Therefore each vortex moves on a cone around the 
momentum vector $. 



We recall that we defined such that el — <1>/||<1>||. From Proposition 3.4, 
it follows that ^ = 0iel for all i = 1, . . . ,2N, hence we obtain the following 
formula for the angular velocity. 

Proposition 3.5 Let Xg = {xi, . . . , X2n) be a relative equilibrium. The angular 
velocity ^ satisfies: 

1 s—^ sin 6i cos 6j — sin 9j cos 6i cos{(j)i — (j)j) 



, ^ 1 — cos 9i cos 9i — sin 9i sin 0, cos(6i — (/>, ) ^ 

for all i = 1, . . . ,2N , where Xi, 9i and are respectively the vorticity, the 
co-latitude and the longitude of Xi. 

3.2.1 Large symmetry relative equilibria 



We first give a theorem which is the analogue of Theorem 3T for relative equi- 
libria. 

Theorem 3.6 Let K be a subgroup of G, x € Fix(J5r) C V and fx = If x 

is an isolated point in Fix(iir) n then x is a relative equilibrium. If in 

addition K is formed of symplectic elements, then x is a fixed equilibrium. 



The proof is an application of the Principle of Symmetric Criticality [P79 and 
the fact that a relative equilibrium is a critical point of the restriction of H 
to a level set The last statement follows from the conservation of the 

momentum map, and the fact that Fyx.{K) is invariant under the flow if K is 
symplectic. 
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Definition 3.7 Let x E V , K its symmetry group in G and ^ — A 
relative equilibrium x is said to he a large symmetry relative equilibrium if x is 
isolated in Fix(X) n <I)~^(/x). 

Large symmetry relative equilibria are relative equilibria with maximal isotropy 
subgroup at a fixed momentum value. Note again that a large symmetry rela- 
tive equilibrium depends only on the phase space, the momentum map and on 
the symmetries of the Hamiltonian: a large symmetry relative equilibrium is a 
critical point for all maps / : <l>^^(/i) M which are G-invariant. 

First, we focus on relative equilibria which lie in a great circle, since their 
dynamics should differ most from planar dynamics. 

Proposition 3.8 For 2N > 8, every large symmetry relative equilibrium which 
lies in a great circle is a fixed equilibrium. 

Proof. To prove this proposition, we could classify all the possibilities by listing 
finite subgroups of 0(3) and examine the dimension of the fixed point set in 
each case. Here is a proof using the properties of the angular velocity. 

Fix a great circle, say {y = 0}. Let x be a relative equilibrium which lies in 
this great circle, and let ^ be its angular velocity. Its symmetry group Gx can 
be written {K, Sy : y i-^ —y) where K is a, subgroup of G. Consider now the 
set of the (+)vortices a;+ = [xi, . . . ,xn) and {K^,Sy) the symmetry group of 
x+ where is a subgroup of 0(3) x 5jv- Similarly, {K^, Sy) is the symmetry 
group of X- = {x'l, . . . ,x'j^). Let tt be the canonical projection t: : K ^ 0(2). 
Without loss of generality, assume that tt{Kj^) C tt^K-). Since Gx — {K,Sy) 
and K+ C K, 

^eFix{{K,Sy),R^) C Fix{(K+,Sy),R^) (Prop. 

Suppose D2 C n{K+), then Fix{{K+, Sy) ,R^) C Fix{{D2, Sy) = 0. Now 
Fix((D2, Sj,), R"^) is the intersection of three pairwise orthogonal planes, hence 
equal to 0. Thus Fix{{K, Sy),M.^) C Fix((i4r+, Sy),R^) = 0, ^ is zero and every 
relative equilibrium is then a fixed equilibrium. 

Cases TT{K-f-) = and 7r(_fir+) = Z2 remain to be treated: in these two cases 
dimFix(G2,) > - 2. Assume now that 2N > 8, then dimFix(G:i) > 2. 
Since a; is a great circle configuration, Fix(Ga;) is a closed submanifold of 
T^^ the 2A^-torus. Hence dimFix(Ga,) n = dimFix(G:,) - 1 and then 

dimFix(G2;) H > 1: large symmetry relative equilibria can not exist. □ 

For 2N < 6, we find easily that large symmetry relative equilibria which 
lie in a great circle are: a (-l-)ring formed of 2 vortices and 2 (— )vortices at 
the poles (denoted G2«(i?, 2p)), two regular 2-rings of opposite vorticities and 
on opposite latitudes with possibly two polar vortices of opposite vorticities, 
and the equatorial (±)rings formed of 4 and 6 vortices (the latters are fixed 
equilibria, see Prop. |3.2| ). 

The following proposition gives a necessary condition to have a relative equi- 
librium formed of two regular A'^-rings. 
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Proposition 3.9 (Lim-Montaldi-Roberts |LMR01|) Consider a relative equi- 
librium formed of two regular N-rings, each of identical vortices with non-zero 
vorticities, together with fcp = 0, 1 or 2 poles, and such that ^ 7^ or /i ^ 0. 
Then: 

- either the two rings are in phase 

- or the two rings are out of phase with an offset equal to . 

We denote these configurations respectively by CNy{2R, kpp) and CNy{R, R' , kpp). 

Explanation about notation Notations of relative equilibria are inspired by 
the Lim-Montaldi-Roberts orbit type notation: that is a relative equilibrium 
Xe = {xi, . . . , X2n) is denoted T{k^Rg,kR,k'R',kpp) where F — iriGx^) (ti" is 
the cartesian projection tt : G 0(3)), R^ is an equatorial (±)ring of 2N 
vortices, i? is a latitudinal A^-ring of identical vortices, R' is a latitudinal Ai- 
ring of identical vortices dual to i?, p is a polar vortex, and k^, k, k' , kp number 
the occurence of Re, R, R' , p in the configuration Xe (the Schdnflies-Eyring 
notation is used for subgroups of 0(3)). 

Now, we show that some of these configurations are relative equilibria. 

Theorem 3.10 Configurations Cnv{2R, kpp) and CNviR, R' , kpp) with kp — 
or 2 poles of opposite vorticities and such that the rings have opposite latitudes 
and opposite vorticities, are relative equilibria. We denote these configurations 
respectively by DNh{2R, kpp) and DNd{R, R', kpp). 




Proof. For each configuration, the method is the following: one computes 
Fix(if) n where K is the isotropy subgroup of the configuration, and 

one uses Theorem |3.6[ For example, consider a configuration Dj^d{R, R') with 
TV odd. Let I be the integer part of N/2. The isotropy subgroup K of such a 
configuration is generated by the two elements 

(r^ o s^, (1 . ..Ny^, (1' . ..N'Y^^t), {Sx,T2.N ■ ■ ■Ti+ij+2,Ty^N' • ••T;/^(,+2)'). 
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The submanifold Fix(K) is parametrized by di and hence Fix(K) n $ ^{fi) is a 
singleton. □ 



Remark. Let N be even and K be the symmetry group of two ahgned semi- 
regular iV-gons of opposite vorticities and on opposite latitudes. K is isomorphic 
to DNh: Fix{K) ~ {(6'i,02) e X S^} and $|Fix(K)(6'i, </'2) = 2iVcos6'i (make 
a picture!). For given fi, —H restricted to V\yi{K)r\<^~^ {n) goes to infinity when 
4)2 goes to (pi or (j)^ which are fixed. It follows that H must have a maximum 
for a certain value of 02- In fact, H is maximal for (/)2 = {4>i — 03)/2 and this 
corresponds to the D]^h {2K) con figuration. This idea can be generalized and 
leads to minimal strata [LMR01|. 



The theorem above does not assume that the polar and ring vortex strengths 
coincide, and in fact it is interesting to be able to vary their ratio. We therefore 
fix the ring vorticities to be ±1 as usual, and let A„ and As be the strengths of the 
North pole and South pole respectively, with A„ — —As. The following proposi- 
tion gives formulae for the angular velocity of relative equilibria _Djv/i(2i?, kpp) 
and DNdiR,R',kpP)- 

Proposition 3.11 The angular velocity of DNhi'2R,kpp) relative equilibria is: 



^ = cos ( 



A^- 1 



N 



l + cosf (j-l) 



(1 + cosf (j-1)) 



The angular velocity of D d{R, R' t kpp) relative equilibria is: 
iV-1 ^ H-cos(f 0--l) + ^) 



^ = cos ( 



N 

^2-sin2 (l + cos(f 0--l) + f)) 



In both cases 6 is the co-latitude of the (+)ring. 

Proof. Let Xe be a Df^[h{2R, kpp) or D^di^^R, kpp) relative equilibrium, and 9 
be the co-latitude of the (-l-)ring in x^. We get from Prop. |3.5| that: 

1 . sin 01 cos 9j — sin 6j cos 9i cos((/)i — 4>j) 



2N+kp 

sin 01 ^ ' 



1 — COS 01 cosdj — sin 01 sm9j cos(0i — ^j) 



Since 



9, A, 



expression becomes: 

N 



^ = COS 

which leads to 
^ — cos 9 



El — cos dj 



1 and Xjy+j = — 1 for all j ~ 1, 

— 1 — cos (j>N+j 

1 — cos 4>j ) 1 + cos-^ 9 + sin^ 9 cos 0Ar+j 



N 



. , iV, this 



sin2 



A^- 1 



sm 



JV 

E 



1 -I- cos (j)N+j 



^ 2 — sin (1 -I- cos07v+j) 



+ 



hp X„ 



cos f sm 
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Figure 3: Angular velocity of Di^ih{2R, kpp). 




Figure 4: Angular velocity of DNd{2R, kpp). 



The formula for ^ follows from (j^N+j = if (j - 1) (resp. ^{j - 1) + f^) for a 
D]\[h{2R,kpp) (resp. Dj\jd{2R,kpp)) relative equilibrium. □ 

For An > 0, the angular velocity of Dnh goes to infinity as the two rings 
approach the equator or poles, while the angular velocity oiDMd goes to infinity 
as the two rings approach the poles and goes to kpXn as the two rings approach 
the equator, we find again that an equatorial (±)ring is a fixed equilibrium. For 
A„ < 0, other behaviour appears for small N (see figures 3 and 4 for A„ = — 1). 
Note that cos d ^~j\ is the angular velocity of relative equilibria formed of a 



single ring of N identical vortices |PD93, LMR] 
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3.2.2 Relative equilibria with less symmetry 

In this section, we give relative equilibria with dimFix(i(r) = 2 for systems of 4 
and 6 vortices, in order to find bifurcating branches of relative equilibria. 

First, we consider configurations of 6 vortices formed of a (+)ring and a 
(— )ring i?_ together with two poles with vorticities A„, As such that As = — A„. 
That is, we consider points of V with symmetry group K — ((r^, ri_2, '''3^4)) 
where "1,2" label the (+)vortices, "3,4" the (— )vortices, and "5,6" the poles. 
Let 9+ be the co- latitude of R+ and 6'_ that of i?_ . The reader can check that 
we can not apply Theorem since Fix{K) n has no isolated points. 

However a critical point of H in Fix(i<')n<I>~^(/i) is a relative equilibrium by the 
Principle of Symmetric Criticality. We therefore use the Lagrange multipliers 
Theorem to find critical points of -ff|Fix(i<') on the manifold *I'|pix(i^ ) (m) order 
to determine relative equilibria. 

For the configurations in question, the Hamiltonian on Fix(_ftr) is deduced 



from (2.1): 

1 2A„ 



H ^In 



((1 — xy)'^ — (1 — .t2)(1 — y^) cos^ a)^ 



{l-x){l + y) 



2{l + x){l-y) 



where x — cos 9+ and y — cos 9^ and a is the offset between R and R' defined 
modulo TT. The momentum map on Fix{K) is given by $ = 2[x — y + A„). 
Applying the Lagrange multipliers Theorem we find that there are two families 
of relative equilibria with isotropy containing K: 

- The first family satisfies a — and 

{y + x)[xy-l + 2A„(y - x) ] = 0. 

The solution y = ~x correspond to the D2h{2R,2p) relative equilibria. For 
|An| > ^, we have a branch of relative equilibria y = (2A„a; + l)/{x + 2A„) 
denoted C2v(2R,2p) and invariant under the symmetry {x,y) )■ {—y,—x) (see 
figure 5). 

- The second family satifies a = f- and 

{y + x)[ x^y^ - 22/2 _ 2x'' + 2xy+l- 2A„(1 - xy){x -y)]^0. 

The branch y = —x corresponds to the D2d{R, R' , 2p) relative equilibria. The 
quartic curve intersects the square [—1,1]^ for all A„ so provides relative equi- 
libria for any value of A„ . That is 



_ X + Xnjx"^ + 1) ± (1 - x^) v^2 + 
^~ a:2 - 2(1 + A„2;) 

is a branch of relative equilibria denoted C2v{R, R' , 2p). Note that C2v{R, R' , 2p) 
is invariant under (x, y) 1— > (— y, —x) and when A„ = under one more symme- 
try [x, y) ^ (y, x). When A„ goes to infinity, C2v{R, R' , 2p) goes to the solution 
y = X that is a (±)ring with two poles of infinite and opposite strength. See 
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-1 ' ' ' • ' ' ' f 




/c2v{2R,2pl 


\ ' ' ' ' .!5 ' ' ' ' 1 

D2h{2R,2p) 









-t ' ' ' ' -.'.5 ' ' ' ' f 


D2hi2R,2p) 



A„ = +1 A„ = -0.6 

Figure 5: Relative equilibria D2h{2R,2p) and C2v{2R,2p). 



figure 6. 

Second, we look for relative equilibria formed of two (+) and two (— ) vor- 
tices with symmetry group K = {(sz '■ z i-^ —z,ti ^4X2 ^3)) where "1,2" label 
the (+)vortices and "3,4" the (— )ones. So we have two longitudinal (±)rings 
Rm,R'm- We proceed as in the previous case: the Hamiltonian on Fix(i4r) is 
given by: 

I 4xy(l — ^(1 — x2)(l — cos a — xy) I 

where x = cos 9i and y = cos 63 and a = ^3 — is the offset between R^ and 
R!^ defined modulo 27r. The momentum map on Fix(ii') is given by $ = 2{x—y). 
We find that there are two families of relative equilibria with isotropy containing 
K: 

- The first family satisfies a = and 



iy+x)[ 2{yx^+xy^-x'^-y'^-xy+l) + {x'^+y^+2xy-2)^yil - x^){l - y^) ] = 0. 

But this does not provide relative equilibria since we have vortex collapse with 
y = —x and the second term does not vanish for {x, y) G [—1, 1]^. 
- The second family satisfies a = w and 

{y+x)[ 2{yx^+xy^-x'^-y^-xy+l)-{x^+y'^+2xy-2)^/{l - x^){l -y'^)]= 0. 

The branch y = —x corresponds to the D2h{2R) relative equilibria. The second 
term corresponds to a branch of relative equilibria denoted C2v{Rm, R'm) which 
is invariant under symmetries {x^y) ^ (—y,—x) and {x,y) 1-^ {y,x). When 
X = y this branch coincides with the Dih{Re) fixed equilibria. See figure 7. 
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C2v{R, R') 




K = 



A„ = 0.6 



Figure 6: Relative equilibria D2d{R, R' , 2p) and C2v{R, R' , 2p). 



Remarks. 

• The equations for relative equilibria Xe with dim Gx^ = n represent the inter- 
section of (n — 1) hypersurfaces of M". 

• One ca n also prove exi stence of relative equilibria using works of Montaldi and 
Roberts [ Mo97 , MoR9£ ] : let Xe be a non-degenerate relative equilibrium (that 
is a non-degenerate critical point) such that ^{xe) = 0, and x an arrangement 
with momentum /i near such that Gx C Gx^- Let be the corresponding 
coadjoint orbit, that is Of^ is the sphere of centre O and radius |/Lt|. One has 
the following: if x is an isolated point of Fix(Ga;, O^), then a: is a relative 
equilibrium. It is an easy exercise to use this result in order to show that there 
exists a relative equilibrium of type C2v{Rrm R'm) near the fixed equilibrium 
D4h{Re) (see figure 7). 



4 Stability of equilibria 

In order to determine stability of equilibria, we use the energy momentum 



method (see Marsden [Ma92| and references therein for a general description). 
Let Xe be a relative equilibrium, /i — <&(xe) G K'^ and f its angular velocity 
(if Xe is a fixed equilibrium, ^ = 0). The method consists of determining the 
symplectic slice G ~ Tx^{SO{3)f_i ■ Xe)^ n Kcr£>$(xe) (so G is transversal to 
SO{3)i_i ■ Xe) where 

50(3)^ = {g e 50(3) | CoAdg -fi = n} 
and then examine wether (P H^\g(xe) is definite or not where 
H^{x)=H{x) + mx)~^i,0- 
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Figure 7: Branches of relative equilibria Z?2/i(2-R) and C2v{Rm, R'm)- 



Note that if d'^H^\g{xe) is definite, that is d^H^\g{Xe) is positive definite or neg- 
ative definite, then Xe is an extremum of the restriction of iJ^ to the symplectic 
slice. 

Let -RT be a subgroup of G. Recall that a relative equilibrium Xe is said to 
be Lyapunov stable modulo K, if for all if-invariant open neighborhoods V of 
K ■ Xe there is an open neighborhood U C V oi Xe which is invariant under 
the Hamiltonian evolution. Papers of Patrick [ Pa92| | and more recently Lerman 
and Singer [1LS98| (see also |OR99| ) refined the energy momentum theorem 
given in | Ma92| . The assumptions of Patrick's Theorem are satisfied since 50(3) 
is compact, we have then the following: 

If d^ H^\g{xe) is definite, then x^ is Lyapunov stable modulo S0{3)fj,. 



For ^ ^ 0, SO{3)fj, is the set of rotations with axis {fi) and then isomorphic 
to SO(2), while 5*0(3)^ = 5*0(3) for /z = 0. Lyapunov stability modulo 50(3)^ 
of a relative equilibrium with non zero angular velocity corresponds to the or- 
dinary stability of the corresponding periodic orbit if 50(3)^ = 50(2), that is 
if /i ^ 0. In the remaining of the section, we will omit modulo SO(3)^. 



4.1 Stability of DNd{R, R', kpp) and DNh{2R, kpp) {kp = 0, 2) 

The most general results of this section are linear instability results: relative 
equihbria of type Dffd{R, R'), DNh{2R), and DNh{2R, 2p) are linearly unstable 
provided N > i, N > 7, and > 9 respectively. The important fact is that 
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linear instability implies nonlinear instability. For A'^ < 6 (resp. N < 8), 
relative equilibria of type D^hi'^R) (resp. D^hi'^R^'^p)) are Lyapunov stable 
for a certain range of co-latitude. First, one needs to describe the symplectic 
slice. 

Notation Let Xe be a relative equilibrium of type D^diR^ R' i kpp) or Djvft,(2i?, kpp) 
with kp = 0, 2, and be the co-latitude of the (-l-)ring of Xg- 

Let j = I, . . . , N number the vortices in cyclic order in the two rings. Write 
tangent vectors {d0j,d4>j) at the j*'' (+)vortex, {d0j,d(l)'j) at the j*'' (— )vortex, 
{6xn,Syn) at the North pole, and {6xs,Sys) at the South pole. 

Let q be an integer, < q < N — 1, and i = \/— T. Define the following 
vectors of T^^V ® iT^^V: 

N 
N 

N 

3 = 1 

N 

where (po is the offset between the two rings, that is </?o = (resp. j^) for a D^h 
(resp. DNd) relative equilibrium. Note that (3q^g, fiqfi', and (3q^^' vanish for 
q = and N/2 (for A'' even) while others do not vanish. 

Define then 5xp = Sxn + Sxs, Syp = dyn + Sys, 6Xp = Sxn — 5xs, and 
^y'p = ^Dn - Sys- 

4.1.1 Description of the symplectic slice 

The symplectic slice G is T^^ {SO{3)n ■ Xg)'^ n Ker D^{xe)- In order to compute 
a basis of Q, one needs expressions of the differential of the momentum map and 
the tangent space to the group orbit at the relative equilibria. 

Proposition 4.1 Let x,, be a DNd{R, R' , kpp) or DNii{2R,kpp) relative equi- 
librium. The differential of the momentum map at x^ is given by: 

D^{xe){66j ,66 J, 54>j ,5(j)'j, 6xn , Sxg , 6yn, 5ys) = 

cos{eo){ai^g + ipi^g) + ism{eo){a[ ,f, + i(3[ ,i,) + Xn{6x'p + i6y'p) ®-sm{eo)a'Q g 

where the direct sum corresponds to the Cnv -invariant decomposition of so (3)* 
as a direct sum of a plane and the line Fix(C„„,5o(3)*). 
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Proof. The expression of the momentum map is 
where 

$r = ^f=i (^i^ '^os <j)j — sin 6j cos (f)'^ , sin 9j sin (fij — sin 9j sin i/)^ , cos 9j — cos 6*^ ) 

Then identifying 5o(3)* with C ® (cz), we get $ = $c © '^z where 

*c - Ef=i (sin Oje^'f'^ - sin ^^-e'*^ ) +A„(x„ - - y,)) 

= Ejli COS0, - cosf?; + A„(Va;2 +yl + + y2). 

Differentiation of this expression at a D^diR, R' , kpp) or DNhi'2R, kpp) config- 
uration gives the resuh. □ 



Proposition 4.2 Let Xe be a DNdiR, R' , kpp) or DNh{'2.R,kpp) configuration, 
and n = $(a;e). 

If fi 0, then 5o(3)^ • x^ is generated by the vector ao,0. 
If fi = 0, then 5o(3)^ • x^ is generated by the three vectors: 

ao,0, Pi,e + cos 9osin9o a[ ,p + Sy'p, ai^g - cos9asm9o + dx'p. 

Proof. Let Xe be a Dnii{R, R' , kpp) or Dffh{2R, kpp) arrangement. The tangent 
space to the group orbit at x^ is: 

5o(3) • a;e = {(77 X xi, . . . , ?7 X X2N+k^) I 77 e M^} 

Hence so(3) • Xe is generated by the three foUowing vectors of T^^V: 

{cx X xi, . . . , X X2N+kp), (cyXXi, . . . ,eyX X2N+kp), (e^ x xi, . . . , x X2N+kp)- 

If ^ = <i>(a;e) 7^ 0, 5*0(3)^ is the set of rotations about the z-axis then 5o(3)^ • Xg 
is generated only by the last one x cci, . . . , Cz x X2N+kp)- The expressions of 
these vectors in spherical coordinates give the desired results. □ 

If 7^ 0, it follows that the dimension of the symplectic slice G is 4:N — 4 
where N is the number of (+)vortices, while if /i = it is of dimension 4N — 6. 



4.1.2 Isotypic decomposition of the symplectic slice 

This second step consists of obtaining a block diagonalization of d^II(^\g{xe) in 
order to compute more easily eigenvalues and then examine definiteness. This 
is done by using the following method of group representation theory (see also 
[ pSS8^ and p^OOC ). 

Let F be a vector space and G be a compact Lie group. Recall that a 
representation W G V oi G is said to be an irreducible representation if W 
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has no proper G-invariant subspaces. There are a finite number of distinct 
irreducible representations of G in V, say Ui,...,Ui. Let 14 be the sum of 
all irreducible representations W C V such that W and Uk are isomorphic 
representations. Then V = Vi (B ■ ■ ■ (B Vi. This decomposition of V is unique 



and is called the isotypic decomposition [3e78|. 

Irreducible representations of groups D^vd and D^h are either of dimen- 
sion 1 or 2. The following two propositions give irreducible representations of 
dimension 1 and 2 respectively. 

Proposition 4.3 Let x^, be a Diq^{R, R! ,kpp) or Dpf}i{2R, kpp) configuration 
and kp — 0,2. There are four irreducible representations of dimension 1 of G^^ 
on T.J.J-' '■ 

(ao,e)K, (ao,0)K, (ao,e)R, (aa0>R • 
Of these, only (ao,e)E and (aj, lie in the symplectic slice. 

Proof. We give only the proof for DNd{R, R' , kpp) arrangements as the one 
for D7v/t(2i?, kpp) is quite similar. The group D^d is generated by two elements 
a,b such that a^ = 1, b'^^ = 1. Let w E T^^V such that a - w — Xw and 
b ■ w — fiw with A,/z e M. Hence = 1, ^u^^ = 1 and A — ±1, ^ = ±1. 
For 1 < j < iV, &^ • {S9i,S(j)i) = {59^(^i),54'c,{i)) where a is a A^-cycle of Sn- 
Thus ((56'o.(i), <5(/>o-(i)) — l^^{^9i,54)i) = {S9i,S(pi). It follows that S9i = S9j and 
S(j)i = S(j)j for all i, j = 1, . . . , A^. The same argument holds also for 69^, 5(j)[. We 
have b- {69i,S(f>i) = {-S9[,S(l)[) = ^i{S9l,S(t>l) then S9[ = ±69i,6<p[ = ±(5(/)i. 
Assume now that 66 1 ^ and 6(pi ^ 0. We obtain a contradiction since 
a ■ {S9i,S(f>i) — {d9i,—6(l)i) = X{59i,6cj)i) and A = ±1. To end, we have to 
deal with polar tangent vectors: a ■ {Sx, 6y, Sx' , Sy') = (Sx' , —Sy' , Sx, —Sy) so 
6x = Sx' , Sy = —6y'. The action of b leads to Sx = 6x' — 6y = Sy' = 0. It 
follows that the irreducible representations of dimension 1 are the ones listed in 
the proposition and only these ones. 

Vectors a^^^ and «[, g do not lie in Q since _D$(xe) -ctQ g — (0, 0, — 2A^ sin((?o)) 
and G so(3)^ • Xe. □ 



Proposition 4.4 Let q be an integer. Assume first kp ~ 0. The following M- 
spaces are irreducible representations of dimension 2 of Gx^ onT^^V: 

* {aq,e,/3q,B), {aq,^,f3q,^), {a'g^g,(3'g g), {a'q^^, P'q^^) , 1 < g < if N is odd 

and Gx, = DNd or Dnh 



* {aq,e,Pq,0), (a,,0,/3g,0), (a'g.e^f^q.e)^ K.^'/^^)' 1 < 9 < A' - 1, 9 7^ f , 

* {ag^e,Pq,e), {aq^ct>,(iq,.j>), {a'^,e^ P'q,e) ^ i^tq,^^ P'q,,p) ' 1 < 9 < ^^1. 9 7^ f, 
(ajv a'jv g), (aiv a'_jv ,) if is N even and G^, = D^h- 
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If kp = 2, then there are two supplementary irreducible representations of di- 
mension 2: {Sxp,Syp) and {6x'p,6yp) . 

The representations (a'l ^, P'l ^) and {Sx'p, Sy'p) do not lie in the sym- 

plectic slice, while the others do. 

Proof. Denote V = Tx^V and Vc — T^^V ® iTj-^V. Irreducible representations 
of dimension 2 are G^;^ -invariant R-planes of V. We can identify M-subspaces of 
dimension < 2 of 1^ with C-lines of Vc by (x, ?;)r i-^ {z = x-\-iy)c where x,y . 
It follows that (a;o:yo)R is a G^;^ -invariant M-plane of V iff xq and yo are not 
colinear and Gx^ acts like the transformation group {{z i— > \z, \z S C}, z i— > z) 
on (zo)c where zq ^ xo + iyo € Vc, that is (2;o)c U (zo)c is a G^^ -invariant set 
of Vc. 

We have Gx^ — i'wd or Djyh- Isometrics in question act on V like: 



{Sx,5y) = {Sx,-Sy) 
{Sx,Sy) = {6x,6y) 

{Sx, Sy) = (cos a Sx — sin a Sy, sin a 6x + cos a 6y). 



Taking account of permutation symmetries, we obtain that the spaces listed in 
the proposition satisfy the above conditions. These representations of dimen- 
sion 2 are irreducible because they do not contain representations of dimension 
1 which are listed in Prop. O. We then use Section [4.1.l| to examine wether 
these representations lie in the symplectic slice. □ 



The following theorem lists the appropriate decompositions of the symplectic 
slice. 

Theorem 4.5 Let N > 3 and Xe be of type D]^d{R, R' , kpP) or Dj^hi'^R, kpp) 
such that /i = '^{xe) is non-zero. Let I be an integer and denote 

Bp = {aa^0,aQ^^,l3i^^,a[g,Sxp,ai^^,l3[g,Syp) 
= (sin 6*0 aifi + cos 6*0 /Jj^^, sin 6*0 Pi^ -f cos 6*0 a\ ,^) 

Bp = ® {ai^e - y cos 6*0 - y cos 6*0 Sy'^) 

Bi = {{aq.e, P'q.^, Pq.e, a'q.4„ a'^.e, Pq..4>^ P'q.e^ o-q:4' I 2 < g < /})• 

Let kp = 0. With the following decompositions of the symplectic slice Q, the 
Hessian d'^Hi\Q{xe) diagonalizes in 2x2 blocks with four 1x1 blocks: 



-kB° ®B^ ® Bn^ for Xe = DNdiR, R') or DNhi'^R), and N odd 
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★ ® ® Bjv_;^ ® (aV fl, /Sjv j,, /3'jv „, ajv >) for Xe = DNd{R, R') and N even 

★ i?° ® i?^ ® Bn__-^ © (o^^.ej Q^'n "^-f ,0' o^'ji 0) •^e ~ DNh{2R) and N even. 

Let kp = 2. FFif/i the following decompositions of Q , the Hessian d^H^\g(xe) 
diagonalizes m 2 x 2 blocks with two 3x3 blocks and two 1x1 blocks: 

-kBl®Bl® Bn_i for Xe = DNd{R, R' , 2p) or DNhi^R, 2p), and N odd 

★ Bl® Bl® Bn_^® for Xe = DNd{R,R',2p) and N 
even 

-k B° (B Bl (S Bn__^ ® (ajv g, a'„ ajv . , a'jv ,) /or Xe = DNh{2R, 2p) and N 
even. 

Wzi/i these decompositions, the matrix L = i7g~^d^7?j|g(xe) of the linearized 
system in the symplectic slice (where fi^ is the matrix of uj\g) diagonalizes in 
4x4 blocks with two 2x2 blocks if kp = 0, and m 4 x 4 blocks with one 6x6 
block and one 2x2 block if kp = 2. 

The case N — 2 will be treated case by case later. 
Proof. There are no difficulties to check that vectors listed form a basis of 



the symplectic slice using Prop. 4.1, 4.2, 4.3, f4.4| and keeping in mind that 



dim 5 = 47V - 4 (resp. 47V) if fcp = (resp. kp^2). 

Let Xe be a Dj^ci[R,R',kpp) arrangement. We have G^^ — C2N >^ '^2[sx\ 
where is the reflection x ^— > —x. Note that C2N is symplectic while is 
anti-symplectic. 

Assume fcp = 0. First, one examines whether some irreducible represen- 
tations of C2N of the symplectic slice Q are isomorphic. The two irreducible 
representations of C2N of dimension 1 are isomorphic. For irreducible repre- 
sentations of C2N of dimension 2, we obtain couples of isomorphic represen- 
tations: subspaces (a^.e, /3q,0, g, a^^^), g > 2 are split in 
(ag,e, Pq,9,a'^^^) and (a'^^g, Z?,,,/,, $'q,s,aq,4,) which correspond respectively to 
rotations of angle qir/N and —q-K/N. The representation B^ is not isomorphic to 
another one, while a'^ g) and l3[ g) are isomorphic. Hence, we get the 

isotypic decomposition of Q, say G = ®jVj with Vq = (ao,e, afj 0): ^1 = B^ and 
dimVj — 4: {j > 2). Since H is G^j^-invariant, it follows that d^H^{xe) is Gx^- 

invariant and the linear map d'^H(^\g{xe) : ®jVj ©j^* is Gxe'Squivariant. 
Let iy^ : Vj G and pv '■ Q* V* be respectively the canonical inclusion and 
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projection. By Shur's Lemma, 

PVj» o d'^H^\g{xe) oivi =0 all i j 

(representations Vj and V* are isomorphic since Gx^ is compact). Since the 
dynamical system is C2jv-equivariant, this result holds also for the linear map 
L = ^t^^ (P H^\g{xe): L diagonalizes in 4 x 4 blocks with two 2x2 blocks. 

The linear map pv o cP H^\g (xe) o iy. is Z2[.S2;]-cquivariant and subspaces Vj 
split in two subspaces of dimension ^ dim V^- (the Z2[sa;]-isotypic decomposition 
of Vj). That is, Vj = [7° ® Uj and 

p^k, o d'^Hclgixe) oiui =0 for all k ^l. 

3 3 

The Hessian diagonalizes therefore in 2 x 2 blocks with four 1x1 blocks. 

Assume kp = 2. With respect to the case kp = 0, the changes are: 

• is isomorphic to the representation {ai,e — ^ cos 6*0 Sx'p, /3i,0— ^ cos 6o ^Up), 
this provides a 4 x 4 block in L and two 2x2 blocks in the Hessian. 

• we have a triplet of isomorphic representations 

{f^i,<l>,<^'i,e)' {<^i,<l>'Pi,e) and {6xp,6yp) 
which provides a 6 x 6 block in L and two 3x3 blocks in the Hessian. 

For D^hi'^R, kpp) arrangements, Gx^ = Cn x Z2[(s2, t)] x Z2[sx]- The proof 
is similar: Cn x Z2 [(sx,r)] is symplectic while Z2[sx] is anti-symplectic. □ 

Remark. The appropriate basis for the block diagonalization does not depend 
on the expression for the Hamiltonian. Let M be the product of 2A'' spheres 
iS^, X G M with symmetry group DncI or Dnh, and : M —^M & Gx-invariant 
map. The Hessian of ^ at a; block diagonalizes like in the case fcp = of the 
previous Theorem. 

4.1.3 Expression of (PH^{xe) 

First, consider the case kp = (no polar vortices). Recall that the Hamiltonian 
is 

H = XiXj ln( 2(1 — cos 9i cos 9j — sin 9i sin 9j cos(0i — 4>j)) ) 

i<j 

and = H + £, J2^^i Xi cos 6i . 

Relative equilibria DNd{R,R') and DNh{'^R) are characterized by coordinates: 

r 9,=9oiov j = l,...,N 
\ 9'j=n-eo{ovj = l,...,N 
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and < 



Ml_llto,j = l,...,N 

f ^ifXe=Z)A..(2i?) 



Note that we assume 9o e]0, ^] without loss of generality since the relative 
equilibria in question have a "North-South" symmetry. Prom now on, we denote 
u = cos 6*0, (fij = ^""^^^^^ and ip'^ = ipj (rcsp. (fj + jf) if Xg is a Djvfc (resp. £>jvd) 
relative equilibrium, in order to lighten formula. 

Some tedious computations show that second derivarives of at 



1 v^A/ 
1-m2 Z^j=2 1 



cos V'j 



J— 2 1 — COSC/2j 

EM — 2«^ + (l — M^) cos <p^ — (1 — cos^ cp^. 
j = l [l + M^-COS^;^ (l-tl2)]2 4^ 



dBidej ~ d0'.ae'- ~ (l-«2)(l-cosvj_jfi) 

aSi^e^ ~ [l+u2-COS<p'j_;^i (1-m2)]2 



8e,90;. - a0i90^. - ^"Vl « [l+„2-cosip^_^, (1-M^)]2 
^H^ _ a^gg _ Y^iV -1 j_ _ 2\^N l-u^-cos¥>^ (iW) 



Pi 

^Hf _ d^Hf 1 



i=2 1-cos ifj + ~ j Z^j=l [l+U^-cosv' (l-«2)]2 



d4>id(pj d<p[d4>'^ 1— cos c^j— frfi ' ^ 7^ J 

''j'-i+i 



k d4>id^'- " )[l+u2-cos¥>' (l-u2)]2 



where the sum in the subscripts is modulo TV. 

Consider now the second case: kp = 2 (with polar vortices). If we permit 
the co-latitude of the relative equilibria to lie in ]0,7r[ instead of ]0, |], we 
can choose A„ = -|-1 and Ag = — 1 without loss of generality. The Hamiltonian 
becomes: 

2JV 

H = AjAj ln(2(l — cos^j cos ^j — sin sin^j cos{(l)i — (j)j)) )+Hr,n+Hr,s+Hn^s 
where 

2N 

Hr,n = '^^i ln( 2(1 - sin Oi cos (j)i x„ - sin di sin y„ - ^Jl-x^- y1 cos 6i) ) 

i=l 
2JV 

ffr,s = - ^ Aj ln( 2(1 - sin 6i cos (t>i Xs - sin 6i sin ys + a/I - a;^ - y| cos 6*4) ) 

i=l 

Hn,s = - ln( 2(1 - XnXs - VnVs + \/\-x\- yl\/l - x^ - yf) ) 
and augmented Hamiltonian: 

/2N \ 



H^=H + i\Y,\,cose. + ^l-xl-yl-^l-xl-yl 
vi=l 
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Note that the angular velocity is different than in the case kp = (see Prop. 



3.5). Second derivatives of i?^ at those given above together with: 

d^H^^ cos d^H^ cos c^J 

d^H^ cos d cosip'^ 



dSidxs l+u Oe'dxs 1-u 

d^H^ sin c^i d sintp^ 

WIdyZ ~ 1-u d9'dy„ ~ l+u 



dOidys l+u de'.dvs 
dtpidxji 1 — u ocp\dxn l+u ' * 

a^ = -^^sin^, = V^sin(^^ 

= -"^^cos^, = "fT^cos(^: 

= i+^cos(p, = -^^^3^ cos 

s^He a^fff a^i/f a^i?; _ 



aa;^ ax^ ay^ ay 

q2 tj ^ 



d4>'dxs 


~ 1-u 




Vl-u^ 


d4'^dy„ 


l+u 




_ Vl-u 


d<t>idys 


1-u 


a%j 




dxjidxs 


dy,idys 



dxndyn dxjidys dxgdyg dxgdy 

4.1.4 Stability results 

One is able now to obtain an expression of (P H^\g{xe) and L = fl^g~^(PH^\g{xe) 
formed of diagonal blocks. We shall see in the proofs of this section that this 
block diagonalization permits us to derive some formulae for the eigenvalues 
of (P H^\g{xe) and L, and thus to conclude about both Lyapunov and linear 
stability. Wc first study the linearization since linear instability is easier to 
prove than Lyapunov stability. Indeed, to prove linear instability, it is sufficient 
that one block is linearly unstable, while to prove Lyapunov stability, one needs 
that the blocks are all positive definite or all negative definite. In the statements 
of this section, it is understood that when a relative equilibrium x is claimed to 
be linearly stable, it means also that the Hessian at x is not definite (thus we can 
not conclude on Lyapunov stability applying the energy momentum method). 
Recall that a linearly stable relative equilibrium can be Lyapunov unstable (see 
[?] for an example). The following theorem gives results for DNd{R, R') relative 
equilibria. 

Theorem 4.6 

• D2diR, R') relative equilibria are Lyapunov stable for 9q G [1.14, 7r/2] and 
linearly unstable on the complement in ]0,7r/2]. 

• D^j^{R,R') relative equilibria are linearly stable on [1.302,1.315] and linearly 
unstable on the complement in ]0, 7r/2] . On the upper bound occurs an Hamilton- 
Hopf bifurcation. 

• For N > A, DNd{R, R') relative equilibria are linearly unstable for all co- 
latitude Oq. 

• Equatorial fixed equilibria D2Nh{Re) are linearly unstable for all N > 3, while 
DihiRe) is Lyapunov stable. 

Proof. Let Xe be a D^diR, R') relative equilibria and > 3. By Theorem 



4.5, the linearized system L = 51 d H^\g{xe) diagonalizcs in 4 x 4 blocks and 
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two 2x2 blocks. All diagonal 4x4 blocks of (PH^\g{xe) and J are respectively 
like (up to a factor): 

/ a c \ 

c 5 

a -c 

V -c b J 

Hence, diagonal 4x4 blocks of L are like (up to a factor) : 

/ c -b \ 

a -c 

-c ~b 

\ a c / 



/ 








-1\ 








1 








-1 








1 1 








/ 



whose eigenvalues are ±i(c + v ab), ±i(c — v ab). These eigenvalues have a non- 
zero real part if and only if a6 < 0. Note that the matrix above is equal to 
cD+aNi+bN2 where D is semi-simple and iVi, N2 are nilpotent. When ab passes 
through zero, linear stability is lost through a Hamiltonian-Hopf bifurcation if 
c 7^ 0. Note also that the corresponding block of the Hessian is definite if and 
only if ab > c^. 

The two 2 X 2 blocks of L correspond to the subspaces {ao^e, ct'o ^) and 
(see Theorem 4^), they are respectively: 





w 



r = (PH^{xe) ■ (ao,e,ao,e) 

) 



where s = d^H^{xe) ■ {a'„ 

w = (fH^{xe) ■ (sin 6*0 ai^e + cos 6*0 /?i_^,sin6'o + cos 6*0 l3[_^). 

The eigenvalues of these blocks are respectively zLiy/rs and ±iw. It follows that 
the subspace is linearly stable and that bifurcation occurs when rs passes 
through zero. 



We first deal with the 2x2 block of L corresponding to (ao,e, aj, ^) 
lengthy computations show that the expressions of r and s are: 



r = 2N 
s = 47V(1 



-{N - l)i±^ + (1 - u4)Ei - 2(1 + m4)I]2 + (1 - u^)j:. 



Some 



where 



and s/N'^ to 



2u 



When N goes to infinity, r/N'^ goes to 



(1 — ) cos ipj 

j^_^2 • One checks that r > for all E [0,7r/2]. One can 
prove that s increases strictly with respect to ^o- Since limits of s at and 7r/2 
are respectively negative and positive, relative equilibria are linearly unstable 
for 60 e]0,tAr[ where is defined by s(iAr) = 0. When N goes to infinity, 
goes to 7r/2. 
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Let N > A. Consider the blocks of L corresponding to the subspaces 
{a'q g, Pq^ip, g,aq^ff,) . FoUowing the notations of the beginning of the proof, 
we have: 



Some lengthy computations show that: 



cos ipj— cos q^j 



a = iV[_(Ar_i) « + 1 v;'^ 

One checks that if q is the integer part of -/V/2, then a6 < for Bq G [t^, 7r/2] 
and for all N > A; L has therefore eigenvalues with non-zero real part, and the 
relative equilibria are linearly unstable. 

For = 3, the matrix L is formed of the two 2x2 blocks with one 4x4 
block, the first 2x2 block is hnearly unstable iff 6*0 G [0, 1.302] and the 4x4 
block is linearly unstable iff 6*0 G [1.315, 7r/2] (the a and h involved are given by 
the above formulae substituting q = I). Since the second block (subspace B^) 
is linearly stable, D^d{R, R') is linearly stable on [1.302, 1.315]. The Hessian is 
not definite on this interval, thus we can not know wether the relative equilibria 
are Lyapunov stable on this interval. 

The D2Nh{Re) fixed equilibria correspond to D^diR, R') relative equilib- 
ria with 6*0 = 7r/2. In this case fi = and 00(3)^ • is generated by ao,0, 



ai,e, and /3i.£/ (Prop. 4.2). The symplectic slice is therefore G = Ker D^{xe) n 
(ao,0j Q^i^e, Hence eigenvalues corresponding to the subspace can be 
ignored when = tt/2. Consider the 4x4 block corresponding to the subspace 

TV cos 

{a[fi,(3i,4„(3[ g,ai^4,): for 6*0 = 7r/2 and q = 1, one has a = ~'2NJ2j=i i^cotU 
and b — N. Since a < for > 3, the D2Nh{Re) fixed equilibria are linearly 
unstable for > 3. 



Consider now the case TV = 2, that is the D2d{R, R') relative equilibria. For 
/X ^ 0, a basis of the symplectic slice is 

("o,e, "0,0, sin 6*0 ai^g ~ cos 6*0 (3[,^, sin 6^0 P[^g - cos 6*0 ai,^). 

The Hessian is a diagonal matrix in this basis. A study of the Hessian and the 
linearized system gives desired results. □ 

Then, consider the D^diR, R' ,2p) relative equilibria. We recall that since 
we assume A„ = +1, wc allow 6q to lie in ]0, tt[. 

Theorem 4.7 
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• D2d{RT R' I'^p) relative equilibria are linearly stable on [2.21,2.31] and linearly 
unstable on the complement in ]0,7r[. 

• D^iiiR-, R' , 2p) relative equilibria are linearly stable on [1.8, 2.05] U [2.25, 7r[ and 
linearly unstable on the complement in ]0,7r[. 

• For N > A, D^diR^R' I'^p) relative equilibria are linearly stable on a small 
interval In and linearly unstable on the complement of this interval in ]0,7r[. 
For TV = 4, . . . , 7, these intervals are: h = [1-75, 1.79], h = [1-73, 1.76], h = 
[1.71, 1.72], It = [1.69, 1.70]; while the length of In is less than 10"^ forN>8 
. Hamiltonian-Hopf bifurcations occur at the bounds of In ■ 

Proof. The proof is similar to that for DNd{R, R')- The components of the 
Hessian are equal to the components of the case kp = except the components 
d^II^{xe) ■ {SOi, 50 i) and of course components involving tangent vectors at the 



poles (see Section 4.1.3). Hence, expressions of s and b do not change, while a 



term — is added to a and a term — ^_ An is added to r. The subspace 

Bp is no longer linearly stable. The expressions of a and h when q = \ do not 
correspond to any eigenvalue formula since a'-^ g and Pi,^ lie in the subspace 
{(3i,^,a[ g,Sxp,ai,jj),f][ g,Syp) which correspond to the 6x6 block. Thus the 
expressions of a and b are useful only for TV > 4. 

For iV > 4, the linear stability is governed by the block corresponding to the 
subspace (a^ g, fiq,ct>, P'q g, oiq,^) where q is the integer part of iV/2 (as in the case 
kp = 0) . Hamiltonian-Hopf bifurcations occur when linear stability is lost (see 
the previous proof). 

The stability results hold also when /i = 0: to obtain the symplectice slice 
when vanishes, one has to remove the subspace 

Q = {Pi,e + cos 6*0 sin 9^ a[^^ + 5y'p, a^fi - cos 6*0 sin 6*0 0^,^ + ^^'j^ 



from the symplectice slice when /i 7^ (Prop. 4.2). For > 4, the subspace 
governing stability has an empty intersection with Q, so the stability of the 
zero-momentum configuration is the stability of this configuration found with 
the generic symplectic slice (where Q is non removed). For N — 2,3, one checks 
by hand that stability results hold also when /i = 0. □ 



Remark. Adding two polar vortices such that the "minus" polar vortex and 
the "plus" ring are in the same hemisphere, tends to stabilize DNd{R, R') rel- 
ative equilibria (above all for = 3): in the sense that the stronger are the 
polar vortices, the larger is the domain of stability. 

The following theorem lists stability results for DNh{2R) relative equilibria. 



Theorem 4.8 

• For N > 7, DNhi'2R) relative equilibria are linearly unstable for all co-latitude 

00- 

• The following relative equilibria are Lyapunov stable in the given range (60 G 
]0,tt/2[) and linearly unstable in the complement: 
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D2h{2R) ]0,0.66] 
Dih{2R) ]0,0.73] 
Deh{2R) ]0,0.45]. 

A pitchfork bifurcation occurs at the upper bound of these intervals. 

• The following relative equilibria are Lyapunov stable on the first range, linearly 

stable on the second range and linearly unstable on the third one: 

D3h{2R) ]0,0.77], ]0.77,0.78], ]0.78,7r/2] 

D5h{2R) ]0,0.67], jo.67,0.68], ]0.68,7r/2]. 

Linear stability is lost through a Hamiltonian-Hopf bifurcation. 

Proof. Assume first that N is odd. The basis for the symplectic slice is the 
same than for D^diR, R') (Theorem 4.5). The eigenvalues of the linearized 
system L on the subspace {ct'n£nPq,<l>T P'g 
where (see Proof of Theorem |4 6| ) : 



9,f ' 



are ±i(c + v a6), ±i{c — v ab) 



b = d^H^ixe 
c — d^H^{xe 



Second, assume that N is even. The blocks are similar to those for N odd 
except the last block: the block of the Hessian d'^H^\g{xe) corresponding to the 
subspace (aw g, a'jv oln_ i, a'jv j) is a diagonal block: 






V 



\ 





b' j 



where < 



a' — d^H^{xe) ■ {piN_ 
a = d'^H^{xe) ■ (a'jv g 
b = d'^H^{xe) ■ {aiv< 



, a N 

2 ' 



b' ^d^H^ixe)-ia'^ 

2 

The associated block of the linearized system L is: 



aiv 



4'' 



/ 





















b 










—a 








V 


-a' 








J 



Eigenvalues of this block are ztiy/ab, ±ivafl/, and bifurcations occur when ab 
or a'b' passes through zero. 

For the two cases N odd and N even, some tedious computations lead to 
the same formulae for a and b: 



- eiVh(iV-l)^ + ^Ef=2 



COS ipj — COS gipj 

1— cos ifj 



EN u —u +(1 — ) cos q'^j-\-{2u —u +1 — (l+w ) cos qi^j ) cos i.pj—{l — u ) cos i 
■\N 1 — cos qtpj 



[l + U-^ — (1 — U-^} COS 

,2^ V^A^ 



\ \ — u —i\+u )cosipj 1 
=2 1-cosipj ' V-^ " J = - COS qipj) [i+„2_(i_„2)(,osyj]^J 

where e = 1 if g 7^ N/2 and e = 2 if g = N/2. One checks that if q is the integer 
part of N/2, then for all 6*0 and all > 8 we have ab < 0, and therefore relative 
equilibria in question are linearly unstable. 
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A study with MAPLE gives results for N = 2, . . . ,7 showing that the blocks 
governing stability are those considered above. The Lyapunov stability is lost 
through a pitchfork bifurcation: D2h(2i?) bifurcates to C2v{Rm, R'm) ^^'^ fo^' 
A^ > 4, DMhC^R) bifurcates to Dn^{AR). For A^ = 2, we have exphcitely de- 



rived the equations of the bifurcating branches in Section 3.2.2. Following the 



method of Section 3.2.2, one can also derive the equations of the bifurcating 



branches for A^ > 4. □ 

Then, we add as usual two polar vortices of vorticity ±1 and obtain the 
following stability results. 

Theorem 4.9 

• For N > 9, DNh{'2R,2p) relative equilibria are linearly unstable for all co- 
latitude 00 . 

• D2h{'2R,2p) relative equilibria are linearly unstable for all Oq. 

The following relative equilibria are Lyapunov stable in the given range and 
linear unstable in the complement in ]0, 7r[; 
D4hi2R,2p) ]0,0.92] 
Deh{2R,2p) ]0,0.83] 
Dshi2R,2p) ]0,0.47]. 

A pitchfork bifurcation occurs at the upper bound of these intervals. 

• The following relative equilibria are Lyapunov stable on the first range, linearly 
stable on the second range and linearly unstable on the third one: 
D3hi2R,2p) ]0,0.83], ]0.83,0.87], ]0.87,7r] 

D5hi'2R,2p) ]0,0.91], ]0.91,0.93], ]0.93,7r] 
Drhl2R,2p) ]0,0.71], ]0.71,0.72], ]0.72,7r]. 

Linear stability is lost through a Hamiltonian-Hopf bifurcation. 



Proof. The proof is a combination of those of Theorems 4.7 and 4.8 



□ 



Remark. The larger is the vorticity of the North pole, the longer is the range of 
stability of D^hi^R, 2p). For example, if A„ = +10 then Dshi'^R, 2p) is stable 
on ]0,1.2]. 

4.2 Other relative equilibria stability results 

In this section, we compute stability of the relative equilibria found in Section 



3.2.2. Relative equilibria in question consist of 4 or 6 vortices. It is not necessary 
to get an isotypic decomposition of the symplectic slice since its dimension 
is either 4 or 8. A study of the eigenvalues of the linearized system in the 
symplectic slice with MAPLE gives the following results: 

• C2v{R, R'), C2v{R, R' , 2p) and C2t,(2i?, 2p) with A„ = ±1 are linearly unstable 
for all co-latitudes. 

• C2v{Rm, R'rn) ^'"^ Lyapunov stable provided the two (-l-)vortices lie in the same 
hemisphere and linearly unstable if this condition is not satisfied. 
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4.3 Bifurcations 



The energy diagrams (see figures 8 and 9) represent the Hamiltonian (the en- 
ergy) with respect to the momentum for the relative equiUbria of systems of 4 
and 6 vortices. Polar vortices are such that A„ = ±1. 

A subcritical pitchfork bifurcation occurs for /i = 1.66, the D2d{R,R') rela- 
tive equilibrium bifurcates to C2v{R, R')- 

A supercritical pitchfork bifurcation occurs for /i = 3.15, the 1)2/1 (2i?) rela- 
tive equilibrium bifurcates to C2v{Rmi R'm)- 

Note that bifurcations are pitchfork since the two bifurcating branches are ex- 
changed by a Z2 symmetry (see Section 3.2.2). 

For the following relative equilibria 1)4/1 (2i?), DehC^R), DshC^R), Deh{'2R,2p), 
D^h{'2R,'ip)i Lyapunov stability is lost through a pitchfork bifurcation: Djsih 
breaks to Divj^. And it seems likely that a second pitchfork bifurcation follows 
immediately: the symmetry breaks to Cjv^ . Following the method of 

Section 3.2.2, one can derive the equations of the bifurcating branches. 



4.4 Conclusion 

The larger is the number of vortices, the more unstable is the relative equilib- 
rium. This agrees with the following well known fact in chemistry: the larger is 
the number of atoms of a molecule, the more unstable is the molecule. 

Relative equilibria DNhi'^R, kpp) are stable near the poles, while relative 
equilibria DNd{2R, kpp) are rarely stable. One can increase the range of stabil- 
ity of these relative equilibria by increasing the vorticity of the polar vortices: 
for DNhi'2R,kpp) one needs to set the "plus" polar vortex and the (-l-)ring in 
the same hemisphere, while for DNd{2R, kpp) one needs to set the "minus" polar 
vortex with the (-l-)ring. 
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